EPJ manuscript No. 

(will be inserted by the editor) 



o 
o 

(N 
c3 



43 

o 

a 

-i— > 
o3 



G 

O 

o : 

(N 
> 

in 



in 
o 
-i— > 

c3 



i 

C 

O 

o 



X 



Intrinsic degree-correlations in the static model of scale- free 
networks 



^ ! J.-S. Lee, K.-I. Goh, B. Kahng a , and D. Kim 



School of Physics and Center for Theoretical Physics, Seoul National University, Seoul 15 1-747, Korea 
fax : +82-2-884-3002 
phone : +82-2-880-1326 

February 2, 2008 

Abstract. We calculate the mean neighboring degree function k nn (k) and the mean clustering function C(k) of vertices 
with degree k as a function of k in finite scale-free random networks through the static model. While both are indepen- 
dent of k when the degree exponent y > 3, they show the crossover behavior for 2 < y < 3 from ^-independent behavior 
for small k to fc-dependent behavior for large k. The fc-dependent behavior is analytically derived. Such a behavior arises 
from the prevention of self-loops and multiple edges between each pair of vertices. The analytic results are confirmed 
by numerical simulations. We also compare our results with those obtained from a growing network model, finding that 
they behave differently from each other. 



PACS. 89. 75. Da Systems obeying scaling laws - 
Self-organized systems 
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1 Introduction 

Structural properties of complex networks have drawn much 
attentions recently 1 1 1213 1 . Degree, the number of edges con- 
nected to a given vertex, is a primary quantity to character- 
ize the network structure. In many real-world networks, de- 
grees are inhomogeneous and their distribution follows a power 
law P c /(k) ~ k~Y. Such networks are called scale-free (SF) net- 
works 1 4 1 . The degree-degree correlation is also important to 
characterize network structure. The correlation between two 
degrees of vertices connected via an edge is measured by the 
mean neighboring degree function k nn (k), which is defined as 
the mean degree of neigboring vertices of vertices with degree 
k 1 5 1 . The correlation among three vertices centered at a vertex i 
is measured through the local clustering coefficient Q, defined 
as Q = 2ei/ki(ki — 1), where e; is the number of connections 
among the kj neighbors of a vertex i. k, is the degree of the 
vertex i. The clustering function C(k) is the averaged one of C, 
over the vertices with degree k 1617 1 . 

While Barabasi and Albert (BA) introduced a model to gen- 
erate SF networks, the model is applied to growing systems 
where the number of vertices increases with time [4|. As an 
extension of the Erdos-Renyi (ER) model of random graph to 
SF networks JS), where the number of vertices in the system 
is fixed, Goh et al. introduced the so-called static model |9|. 
The term 'static' originates from the fact that the number of 
vertices N is fixed. The static model was followed by other 
similar-type models such as the hidden variable model ifTol 
II 11121 . In the static model, each vertex i(i = ,N) has a 



prescribed weight pi summed to 1, which is given as 



(1) 



where the Zipf exponent /j is in the range < /j < 1 . To con- 
struct the network, in each time step, two vertices i and j are 
selected with the probability pi and pj, respectively. If i = j or 
an edge connecting i and j already exists, do nothing, implying 
that self-loops and multiple edges are not allowed, respectively. 
This condition is called the fermionic constraint hereafter. Oth- 
erwise, an edge is added between i and j. This process is re- 
peated NK times. The resulting network is a scale free one with 
the degree exponent given as 191 131 
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Since a pair of vertices is selected with rate Ipipu where the 
factor 2 comes from the two cases of and one may 
think that there is no degree correlation, which is the case we 
can observe when y > 3. However, when 2 < y < 3, due to 
the fermionic constraint, the degree-degree correlation arises 
intrinsically. In this case, the degree correlation occurs for ver- 
tices with large degree, while it is still absent for vertices with 
small degree. In this paper, we investigate such degree corre- 
lations in terms of the functions k nn (k) and C(k) in the static 
model and their crossover behavior in terms of system size N. 

Many SF networks in the real-world and artificial networks 
include degree correlations within them. For example, the mean 
neighboring degree function k nn (k) behaves ~ k~ v with v > 
for the Internet |5| and the protein interaction network |14|, 
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while v < for social networks such as the coauthorship net- 
work. The case with v > (v < 0) is called disassortative (as- 
sortative) mixing |15|. When a network contains hierarchical 
and modular structure within it, it is suggested that the mean 
clustering function C(k) behaves as C(k) ~ k~$ for large k as 
observed in metabolic networks and the hierarchy model |6 
1161 . Occurrence of such degree correlations in real- world net- 
works may be related to their own functional details. For ex- 
ample, the assortativity of the social network arises from the 
social relationship between bosses, while the disassortativity of 
the Internet comes from the network design to allow data pack- 
ets flow efficiently. The three-degree correlation may be related 
to the control system in biological network such as the feed- 
back or feed-forward loop structure 1171 . Such degree correla- 
tions in real world networks appear in the combination of those 
due to the fermionic constraint and their functionalities. On the 
other hand, the static model is frequently used to study vari- 
ous dynamical properties of complex networks. Therefore, the 
knowledge of the intrinsic degree correlations we study here 
would be helpful in understanding the degree-correlation a SF 
network has for functional activity. For the purpose, Catanzaro 
and Pastor-Satorras |18| studied the degree-correlation func- 
tion k nn (k) for the static model, but their study relies on numer- 
ics in the end. Here we present analytic solutions for k nn (k) 
as well as clustering function C(k). We mention that k nn (k) 
for a related model was analyzed by Park and Newman 1 19 1 
while k lm (k) and C(k) for the BA-type growing network mod- 
els are studied by Barrat and Pastor-Satorras 1 20 1 using the rate 
equation approach |21 1. On the other hand, it was desirable to 
introduce uncorrected SF network as a null model to check 
the correctness of analytic solutions in various problems on SF 
networks. For the purpose, Boguna et al. |22| and Catanzaro 
et al. [23 1 introduced a way to construct uncorrelated SF net- 
work by restricting degree of each vertex to be less than the 
cutoff value k c , beyond which the intrinsic correlation arises 
in 2 < y < 3. The cutoff value they used scales as ~ N x l 2 , in- 
dependent of y, which was based on the configuration model 
introduced by Molloy and Reed |24|. Such cutoff is also im- 
plicit in the model introduced by Chung and Lu 1251 . How- 
ever, we show that while the natural cutoff of the static model 
is ~ jyVOr- 1 ), the vertex correlations appear for degrees larger 
than a crossover value, k C[ ^N^ 2 ^^ 1 ^, which is smaller than 
N x l 2 for 2 < y < 3. We mention that k C[ — consul N l l 2 ) as 
y — > 2(y — > 3) so that for y — > 2 all the nodes have nontrivial 
vertex correlations and for y — > 3, there are no correlations. 

In Section 2, we derive the mean neighboring degree func- 
tion k nn (k) and the mean clustering function C(k) analytically. 
Comparisons between the results of our analytic derivations 
and numerical simulations are given in Section 3. Section 4 
summarizes our results. 



2 Analytic Results 

In the static model, the notion of the grand canonical ensem- 
ble is applied (13)-, where the number of edges is a fluctuating 
variable while keeping the SF nature of the degree distribution. 
Each pair of vertices (z, j) is connected independently with the 




Fig. 1. The connection probability fij of an edge has two distinct re- 
gions where fij ~ 1 or ~ 2KNpipj due to the fermionic constraint 
when 2 < y< 3. 

probability fij, given by 

f ij = \-e- 2NK W, (3) 

because the probability that vertices i and j (i ^ j) are not con- 
nected after NK trials is given by (1 - 2piPj) NK ~ e" 2NKpip i . 
That is, if we denote the adjacency matrix element by fly (— 
0,1) then its ensemble average is given by fij; i.e., (ay) = 
fij, (■ ■ ■ ) denoting the grand canonical ensemble average. For 
i = j, fij = because of the prevention of self-loops. Since 
INKpiPj ~ KN 2 »- l /(ijY for finite K, when < p < 1/2, cor- 
responding to the case y > 3, IKNpipj is small in the thermo- 
dynamic limit, therefore, 

fij « 2KNp iP j. (4) 

This is called the bosonic limit. On the other hand, when 1 /2 < 
p < 1, corresponding to the case 2 < y < 3, IKNpiPj may di- 
verge in the thermodynamic limit, therefore, fij is not neces- 
sarily of the form of Eq. @, but it reduces to 

, / 1 when ij<^iN 2 ~K ... 

fij ~ \ i PJ 

{ IKNpipj when ij > N f. 

This is the manifestation of the fermionic constraint, the pre- 
vention of multiple edges. Thus, for 2 < y < 3, fij has two dis- 
tinct regions in the (z'j) plane as shown in Fig.^ 

2.1 Degree and degree distribution: 

The degree of a vertex i is given in terms of the adjacency 
matrix as ki = £ ; fly. For completeness, we present here known 
results for the mean degree (kj) [13|. It is obtained through 
the formula (kj) = Y*jj=ifij which can be evaluated by using its 
integral form as 

u- 1 

rN 1 1 i 1 f aN 1 1 — e~ xy 

(k,)=Y,f<J- djfij = -a*N' * , dy — , 

h p JaN~7 y< 

(6) 
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where y = 1 + -, x = aN^ 2 i M and y — aW 2 j p with 



a=^2K{\-nf, (7) 
The integral in Eq.0 denoted as I(x) is evaluated as 

t dtf~i i — e -*y 

I(x) = __ x dy- 



laN 2 



a 2 --lN~ 



y 



7 



Y-2 
y-2 



-(l-N> J - i )x for x<l/aN»-?, 
-x + qo {y)x y ~ 1 for x > 1 /aW- 2 . 



(8) 



with ^o(y) = Jo ' dr( \ —e r — r)/r*, which is a negative con- 
stant. Therefore, we obtain that 



(ki)^2K(l-fi)(-r) +SL, 
where A is the correction, of which the leading term is 



(9) 



_ f -2isT(l -rfN 2 *- 1 /? for i>aT'N 2 K 

\(2K(l-ij) 2 y'N 2 ~vq {Y)/(iij) for i<aT-N 2 ~T.. 

(10) 

This is negligible compared with the first term on the right hand 
side of Eq. (|9j in the thermodynamic limit, N — > °°. The average 
degree is then 



2(L) = 



(11) 



where (L) is the mean number of edges in the grand canonical 
ensemble. From Eq.(|9j, one can easily obtain that the degree 
exponent is related to the Zipf exponent /j as y = 1 + 1//J given 
inEq. (0 



2.2 Mean neighboring degree function k nn (k): 

We pay attention to the case 2 < y < 3. To evaluate the mean 
neighboring degree function k nn (k), we first calculate the mean 
neighboring degree in terms of i, i.e., k nn (i) and convert it to 
k n n(k) by using the relation of (k() versus i. To proceed, we use 
the expression, 



k\m(i) 



Ej€im of i^j \ / Y l j,k a ij a jk \ (T,j,k a ij a jk) 



Eyfli. 



(ki) 



(12) 

where the ensemble average is applied to the numerator and 
the denominator separately. Its validity is checked numerically, 
which is shown in Section III. The denominator was already 
derived, and the numerator is evaluated as follows: 



Y, a 'j a jk) = L f<jfjk + Y,fij ~/ d J dkfijfj k +{ki). 



(13) 

where a?- = a i; - is used and the double sum is approximated 
by the double integral. The validity of the transformation from 
the discrete double sum to the double integration is discussed 



in Appendix A where 2 < y < 3. Such an approximation intro- 
duces at most an 0(1) factor on the amplitude of the leading 
order terms for large N, as will be mentioned below. Applying 
the similar method used in Eq. 0, we evaluate the integration 
in Eq. Jl 3i as 



dj / dkfijfjk - 

aA^~2 i_ e -xy t-aN"~7- \ _ e ~yz 



1 dy — it - 

aN 2 y> JaN 2 Z> 



dz- 



(14) 



where x = aN fJ ~ 1 / 2 r fJ is in the range aN~ x l 2 < x < aN^ l l 2 
(see the definition of a in Eq.flTJ). The last part of the integral 
in Eq. (114-i is I(y) defined in Eq. 0. Therefore, we substitute 
the leading term of Eq. (|8} into Eq. J14I and obtain 



dj / dkfijfjk w — xi / , dq . 

1 Jl fi(l ~fJ) JaxN~T- q< 1 



(15) 

in which we change the variable of integration as q = xy. The 
integral in the right hand side of Eq. Jl 5I > is evaluated in three 
parts as 



axN" 2 l- e -9 /■<*> l- e -1 

dq — : — = / dq- 



axN~ 2 qi 1 Jo q y 1 



axN 2 j _ e -q 

^ q<- 1 JaxNK-h"' 1 q1- X 



The first term is denoted as 

91 (Y) = 



1 - e-'i 



q 



1 «9 — ^n— -(16) 



(17) 



which is a positive constant for 2 < y < 3. The second term is, 
since axN~ l l 2 < 1, 



/„ 



axN 3 l-e^ a 3 -TAr-( 3 -Y)/2 



9 



,7-1 



3-y 



■3-7 



(18) 



_1 2 2-I 

The last term is calculated as, when x -C 2 (i^a^N ^), 



dq——r- = \ dq—^ / dq- 



axN^i q y 1 Jo <? T 1 Jo 



91 (Y) 



a 3-Y#(3-Y)0H) 



3-y 



(19) 



and, when x > 1 /aiV 2 (i^a^N f ), 



1-1 



.2-7 



(20) 



Combining all the contributions, when x -C 1/ flA* 2 , the sec- 
ond term on the right hand side of Eq. d!6i becomes the leading 
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order term, while when x » 1 /aN p 2 , the first one does. Thus, 



The numerator is evaluated as 
1 



for i>a^N' 



2- 



dj dkfijfjk'' 



g4jy3g-l 
(1-^-1)' 

?i(y)4(Sr rl+ ' i for /<fli ^ 



5 



E f'jfjkfi 



ki 



(21) 

The second term (&,•) on the right-hand-side of Eq. Jl 3i can be 
neglected compared with Eq. ( 12 li for all range of 2 < y < 3 and 
i. Therefore, 



^nn (0 



- 2-i 

when i>afN '', 



<?i(y)a" 2 JV 3 2 ^ ~n 2 >-'- 1 /^ when i<aT-N 2 



and using Eq. (|9} for A: = (A:,), 



(22) 



\ ^-rN 2 **- 1 when k<N l ^ : 

M*)«i ^ 1 2 _1 2+ l 

[?i(y)2A'(l-|i)/!Ar *!Jt when k>N 1 ^. 

(23) 

?-i 2 2- A 

Here we note that the coefficient of N p when i > at'N '< (or 

when k < N l ~ >J ) is not exact but is in between a 2 / (2/j — 1) and 

2/Lta 2 /(2/j— 1) as explained in Appendix A (see Eq. (l44t ~). In 

terms of the degree exponent y we rewrite Eq. J23l > as 

k f«~/ M3 ^ /( r^ when k>kci > (24) 

where the crossover degree k c \ scales as k cl ~iV(Y-2)/(Y-l). 



2.3 Clustering function C(fc): 

The clustering function C(fc) is the mean of the local cluster- 
ing coefficient C, over the vertices with degree k. To calcu- 
late C(k), we first calculate C, and convert it to C(k) by us- 
ing the relation Eq. (|9}. As we introduced before, C,- is defined 
as C; = 2ej/kj(ki — 1), where <?,■ is the number of connections 
among the A, neighbors. In the grand canonical ensemble, Q is 
calculated as 



Q = 



ki(ki-l)/2 



However, we approximate it as 



Q 



(et) 



(ki(ki — \)/2)' 



(25) 



(26) 



which enables us to calculate it analytically. The validity of 
this approximation is checked numerically in Section III. We 
evaluate the denominator and numerator separately. 
The denominator is evaluated as 



ki{ki-\) v _l v , , 

\ L fijf ik ™2K 2 {\-n) 2 N 2 »/i 2 ». (27) 



^Lfi}fjAi~2h dj Jl '' kfnlikfu 2^ 2 " 



aN 2 



aN 2 



Possible errors involved in using the integral form for the dou- 
ble sum is discussed in Appendix A and will be mentioned be- 
low. The evaluation of the integrals of Ea. J28t is carried out 
depending on the magnitude of x. When x 3> l/aiV _ 2 , we ob- 
tain 

raN"-h. (1 _«-*)(! 
/ 1 dz v 

f 40 (y) (xf~ l -(x+y) y - 1 ) when y < 1 /aW" 2 , 

~ \ 90 (Y) (* y_1 + 3' Y ~ 1 - (x + y) 7 ^) when y > l/aN^. 

(29) 

Thus (e,-) is written as 



1 2- A 



2^2 



9o(y)(® + c), 



(30) 



where $ and C are expressed in the integral forms as 



raN^ 2 1 _ e -xy 
<B = I , dy — 



aN 
1 



dq 



(y-l)(l-e- xe >) 



.7-1 



I <fy 



oW 2/jt 
aN^'l/x 



y-1 1 



(1 



and 



C = -?o(Y) / i 

/a AT 2 



l/a2V* 2 J _ e -xy 



dy- 



y< 



-y 



(31) 



(32) 



Even in the region of x > 1 /aN 11 2 , the leading term is de- 
termined depending on the magnitude of x. When x > 1 > 

l/aN^?, the first term of the integral $ is the most domi- 
nant one compared with the other terms of <B and C, which is 

evaluated as w (y- \)q\{y)x 2 ^~ 2 \ When 1 > x > l/aN p ~?, 
however, the third term is most dominant and evaluated as « 
2(y — l)ln(l jx)x 2 . Therefore, the numerator is evaluated as 



%I ^ L [-qo(i)qimj- \)x 2{ ~l-\ when x > 1, 

A ' 2(y-l)ln(l/x)x 2 , when 1 > x> l/oiV - 2. 

(33) 
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Thus, we get C, in the region of x > 1 /aN M 1 to the leading 
order as 

c _ (N 1 ~?\n(P/N''-z), when N 1 ^^ < j<Ar 2_ ^, 
S N 5-fy-l. 2(2^-1) ^ W hen/<A? 1_ 5;. 



Equivalently, 



C(t) 



?\n(m/k), when N 1 - 1 " <k <N?, 
A7 3 "Ijfc-2(3-Y) ) w henAT2<A:. 



(35) 



1 2— - 

Let us consider the case of x < 1 /oA"^ (i > AT ?). In this 



case, 



off i (i_ e -Ji)(l_ e -«) 



_ i 

aiV 2 



for y > 1/ aN M 2 and is almost negligible for y < 1/ aA^ 1 2 . 
Thus 

2 1 1 

anMavN^v f aN ''~ J , , l—e^ 

qoW a I [( x + 37 )T-l_ 37 7-l] i 



2^ 



g o(y)qMv-V-^ 

2jU 2 v y 



(37) 



_1 2 2-1 

Therefore we get when jc <C 2 (i.e., 3> a^Af f , or 

Q=C(ifc) ^AA^lnAT, (38) 

with A = -^ (Y) fl ^ 2 (Y-2) 2 (3 -y), when k < N 1 ^. In Ap- 
pendix A, the error introduced in Eq. d28l > is estimated and is 

found not to change Eq. ( 1351 . For Eq. ( I38K however, we find 
,i 

C(k) ~ AN * (lnAT + D) with an undetermined constant D. 
Eqs. J35l > and d3 8I > are the main results of this subsection. 



3 Numerical simulations 

We now discuss numerical check of the analytic results derived 
in Section 2. For the purpose, the static model network is gener- 
ated with K = 2 and fu = 2/3 (y = 2.5) and with varying system 
size N. All data below are averaged over 10 4 network config- 
urations. For the case of k m (k), we first check the approx- 
imation, Eq. M2\ . To proceed, we measure L/e nn of /(fc/Ai) 
and £,- 6 nn of i (kj) / (£;) separately in Fig|2l finding that the data 
overlap and the approximation is valid. Next we directly enu- 
merate the function, J^ 1 d j J r dkfijfjk / d jfij + 1 (solid line) 
and compare it with the evaluation (dashed line) within leading 
order, Eq. d23i . The extra term of ' 1' comes from the 2nd term 
(ki) of Eq. Jl 31 . For small i, the two lines seem to be consistent, 
however, for large i, they somewhat deviate in the intermedi- 
ate region of i. However, we confirm that our analytic solution 
is valid within leading order by the finite size scaling plot. In 
Fig|3 we plot k nn (i) for different N = 10 3 , 10 4 , 10 5 and 10 6 
finding that the data collapse into a single curve by the rescal- 
ings of i — > i/N 2 ^ 1 ^ and k nn (i) — > k nn (i)/N 2fJ ~ l . Moreover, we 




10' 



„ -iisnnofi < *i''* : i* > + 

S j is nn of ,<*;>'<*', > O 

Integration 

Asymp. tor small i — 

plateau for large i 



10" 



10' 



10* 



10" 



10° 



10° 



Fig. 2. Plot of fcnn(s) versus i. To check the validity of the approx- 
imation Eq. A3, we plot £j 6nno n(|) (+) and £ ' e " n " f i{kj) - 



(36) ^ijfr'* + 1 (°) f° r N = 10". We can see that the approximation is 



valid. We compare them with ^-jg^fe^ + 1 (solid line). They agree 
with each other for small i, however, it is in disagreement in the plateau 
region as expected in Appendix A. We also plot with the first leading 
term presented in the text in the asymptotic regions with dot-dashed 
line and dashed line, respectively. 



5 




A*^ 10" 



10-' 



10" 

N 



Fig. 3. Plot of size-dependent behavior of k nn (i). Data of different 
network sizes N = 10 3 (D), 10 4 (*), 10 5 (a) and 10 6 (4) collapse into 
a single curve in the scaling plot. Inset: Plot of the difference between 
the leading order analytic expression and the simulation value of i = 1, 
divided by the simulation value. As N increases, the relative difference 
decreases, showing that the analytic solution converges to the numer- 
ical result. 



find as we increase the system size that the numerical simula- 
tion data approaches our analytic solution for small /(inset of 
FigQ. We also check the behavior of k nn (k) numerically. Un- 
der the rescaling of k -» k/N 1 ^ and k nD {k) -> k n n(k)/N 2 ^ 1 , 
the data for different system sizes collapse well, confirming the 
validity of our analytic result. 

Next, the local clustering coefficient function C, is mea- 
sured. We first check the approximation introduced in Eq. fl26> 
in Fig|5] finding that they overlap each other except for large i. 
This discrepancy originates from the fact that the vertices with 
large i are mostly those located at dangling ends with degree 
1 . Thus, the formation of triangles or wedge shapes is rare and 
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^ 10" 



10" 




10' 



10" 



10° 10 



10' 



10° 



Fig. 4. Plot of k nn (k) versus k for different system sizes N = 10 3 (D), 
10 4 (*), 10 5 (a) and 10 6 (^). Data for different system sizes collapse 
in the scaling plot. Solid and dot-dashed lines indicate the analytic 
results of leading order for large and small k, respectively. 




JL5L 



o o 



<eflk i (k r l)/2]> 
<e j >l<k i lt r l)/2> 
Integration 
Asymp. for small i 
plateau for large I 



10° 



10° 



2e, 



Fig. 5. Plot of Cj. We can see that the approximation, {^j- 

(Hkj-\)) ' s va 'i^ f° r * ne large ' limit. The dot-dashed line indicates 
the analytic result. In the plateau region, the discrepancy between the 
analytic and the numerical results decrease as system size N increases. 
The dashed line indicates the analytic results, Eqs. 1341 and I38i . 
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10 10 



10" 

N 



10 10 



10" 10 

.(Y-1)/(6-2y)-, 



10' 



10 J 



10"' <-T- — ^ 

10" 4 10" 3 10' 2 10" 1 

Fig. 6. Size-dependence of local clustering coefficient Q. Data of var- 
ious network sizes N = 10 3 (D), 10 4 (*), 10 5 (a) and 10 6 (» are col- 
lapsed in the rescaling plot. Inset: Plot of the difference between the 
analytic solution within the leading order and the simulation value for 
1 = 1, divided by the simulation value as a function of N. The decreas- 
ing behavior with increasing N indicates that the analytic solution is 
asymptotically valid. 



plateau 
10 3 , 

10= 
10 6 



10- 




Fig. 7. 

10 5 (A 



10 10 10 u 10 10 10 

Plot of C(k) for different system size N = 10 3 (D), 10 4 (*), 
and 10 6 (4). The data are well collapsed in the rescaling plot. 



their numbers fluctuate highly. Next, we also check the valid- 
ity of the approximation from the discrete summation and the 
continuous integration, Ea. J28l l. For small i, the approximation 
is reasonably valid as shown in Fig|3] which can be expected in 
Appendix A. However, for large i in the flat region, the approxi- 
mation shows some discrepancy, but it is likely that the discrep- 
ancy decreases as system size N increases. To check the size- 
dependent behavior of C,-, we plot C, versus i with rescalings of 
Q -> Q/N l -hnN and i -> // (^(^-4+^/(4^-2) ln i/(4p-2) ^ 

for different system sizes N = 10 3 , 10 4 , 10 5 and 10 6 . We find 
that the data collapse reasonably well as shown in Fig|6] And 
we also check the behavior of C(k). By rescaling of C(k) — > 

C{k)/N l ~* \nN and k -> Hn 1/2(2 ~^ /JV 1 / 2 , the data of C(k) 
for different system sizes also collapse into a single curve rea- 
sonably well as shown in Fig0 Thus, our numerical simulation 
results show that, although several approximations are involved 
in deriving the analytic results of section 2, they are valid to the 
leading orders in N as N — > °o. 



4 Conclusions and discussion 

We have studied analytically the mean neighboring degree func- 
tion k nn (k) and the clustering function C(k) in the static model 
for the case of 2 < y < 3 and checked the results by numerical 
simulations. Due to the prevention of self-loop and multiple 
edges, there occur intrinsic degree correlations, which appear 
for 2 < y < 3 in the fc-dependent form of k nn (k) and C(k) for 
large k. Our results are summarized in Table I together with 
those for the case of y > 3. It would be interesting to com- 
pare our results with those obtained in the generalized BA-type 
growth model [201 ■ In this model, k nn (k) ~ Ar(3-Y)/(Y-%-(3-Y) 
when y < 3, ~ InN when y = 3, and ~ Ink when y > 3. On the 
other hand, C{k) ~ aK4-2y)/(y-%-(3-y) for k > 0niV) 1/{3 "^ 
and ~ (lnAOM 4 - 2 ^- 1 ^-^ 3 ^) for k < (lnA?) 1 /(3-Y) whe n 
y < 3, ~ (lnN) 2 /N when y = 3 and ~ N~ l k?~ 3 in the range 
k < Af 1 /^ -3 ) when y > 3. Therefore, it appears that the degree 
correlation functions k nn (k) and C(k) behave differently for the 
cases of the static model and the BA-type growth model. 
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Table 1. Degree and system-size dependence of k nn (i), k m (k), C; and C(k). 



2<y<3 



Range 


knn(i) 


knn{k) 


Q 


C(k) 


i > N 3 -y 

k < Ni- 1 


3-7 

~ A/y-i 


~ AMY" 2 ) \nN 


N 3 -1>i>N^-t 
Nr^ <k<N? 


(3-Y)(y-2) 3-y 


~N 3 -yk-( 3 -y') 






Arl(3-Y) > i 

Ni <k 


~ N i- 1 it- 1 


2r+9y+ll 2(3-T) 

~ A' i y-' 


~ N 5 - 2 ^k- 2( - } -^ 



Y=3 



whole range 



lnTV 



(lnN) 2 /N 



Y> 3 



whole range 



-0(1) 



UN 



This work is supported by the KRF Grant No. R14-2002-059-010000- of the surface term to (k lm ) (i) is - a 2 N 2 ^ 1 . Then Eq. lf22|l for 
in the ABRL program funded by the Korean government MOEHRD. ^ 



- 2-i 

: case i > a^N '' has to be changed as 



A Transformation from discrete summation 
to continuous integration 

In several parts of this paper, we use the transformation from 
the discrete summation to the continuous integration such as 



fi dj fi dkfjjfj^ < < Sfdj Si dkfijfjk t 



(ki) 

This leads to 



(ki) 



a 2 N 2 M -\ _ a 2 N 2 »- 1 



(ki 



(43) 
(44) 



j,k=i 



N rN 

dj / dkF(i,j,k). 



Here we discuss its validity. For a monotone decreasing func 
tion F(x), one has the well known relation: 



(39) to the leading order in N. Thus, the leading order of (k nn )(i) is 

- 2— - 

given only in the form of the bounds when i> a^N '' . 

When Ci is considered, F(i,j,k) = fijfjkfki- The most rel- 
evant terms are Y,kF(h ^->k) and Y,jF(hj> 1) : 



r N „n 

dxF{x)+F(N) < Y F(n) < / dxF(x) +F(l) 



£F(i,l,*)=£F(i,j ) l) 



(40) 



Y^filfikfki « -«o(Y)(Y- l)ah 2 -^x 2 /^ (45) 



When F(i,j,k) is positive, monotonously decreasing and bounded 

in both /' and k, we can apply Eq.Eol twice to obtain the error , . _ t,* t 2-t, tt. tt) , j e t A 

. -J— ' FF J ^' — ' when;>a''V ^ . These are of the same order of magnitude as 

in Eq.(|39) as 



Y,F(i,N,k)+Y,F(i,j,N)-F(i,N,N) 

k j 



<Y,F(i,j,k)- dj dkF(i,j,k) 
< \,k) +T^F(i,j, 1) -F(i, 1, 1). 



(41) 



Thus Ea.J39> is valid when the "surface terms" in Eq. (l41l are 
negligible compared with the "bulk term", Si d j Si dkF(i 7 j,k). 
When we consider k nn (i), F(i,j,k) is given as fijfjk and one of References 
the surface terms that require special attention is 



the bulk term in Eq.(l33> up to the ln/V factor. Other boundary 
terms are smaller in order of magnitude compared with these 

2 2-1 

terms. Thus when i > a^N f , C; is bounded as 

- ? o(Y)a^ 2 (Y~ l)(Y-2) 2 A^ 1_ hn(fl 2 A^ 2 ^ 1 ) < C, 

< -q (y) a y 2 {y-l)(y-2) 2 N l ~» [ln^/V 2 ^ 1 ) + 4,u] . (46) 

The boundary term is important when ln/V is not large enough 
compared with 4/u. 



£F(M,*)=£/a/i*: 

*: k 



-1 



(42) 



with a = ^2K{\ — /j) 2 . It turns out that this surface term is 

2 2-1 

of the same order as the bulk term when i > wN f. Other 
surface terms are, however, negligible. Thus, the contribution 
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